Elastic proton scattering of medium mass nuclei from coupled-cluster theory 
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Using coupled-cluster theory and interactions from chiral effective field theory, we compute overlap 
functions for transfer and scattering of low-energy protons on the target nucleus 40 Ca. Effects of 
three-nucleon forces are included phenomenologically as in-medium two-nucleon interactions. Using 
known asymptotic forms for one-nucleon overlap functions we derive a simple and intuitive way of 
computing scattering observables such as elastic scattering phase shifts and cross sections. As a 
first application and proof-of-principle, we compute phase shifts and differential interaction cross 
sections at energies 9.6 MeV and 12.44 MeV and compare with experimental data. Our computed 
diffraction minima are in good agreement with experiment, while we tend to overestimate the cross 
sections at large scattering angles. 



PACS numbers: 21.10.-k, 24.10.Cn, 24.50.+g 

Introduction - With the advent of accelerators of new 
generation, providing radioactive ions beams, it becomes 
possible to synthesize nuclei far from the valley of stabil- 
ity. Little is known about these nuclei. For example, the 
nuclear interaction is not well understood at large proton- 
to-neutron ratio, and it appears that its effects are radi- 
cally different close to drip-lines compared to the vicinity 
of the valley of stability. The most striking example is the 
redistribution of shell closures, which have been noticed 
to be different from the usual magic numbers present in 
well-bound nuclei PQH]. In order to precisely study nuclei 
close to drip-lines, new theoretical tools must be devel- 
oped, going further from the standard structure methods 
of nuclear analysis, based on the use of standard shell 
model, treating all nuclear states as well bound, and, in 
reaction theory, of optical potentials fitted from exper- 
imental data. Indeed, both nuclear structure and reac- 
tions must be unified, so that both inter nucleon cor- 
relations and scattering degrees of freedom are included 
within the same framework [3] . Microscopic models aim- 
ing at a unification of nuclear structure and reactions 
have made great progress over the last decade. The Res- 
onating Group Method within the No Core Shell Model, 
has successfully described nucleon and deuteron scatter- 
ing and fusion in light nuclei @], the Green's Function 
Monte Carlo method has been used to describe elastic 
scattering on 4 He and in the computation of asymptotic 
normalization coefficients in light nuclei [5], and finally 
the Self Consistent Green's Function (SCGF) method has 
been applied to the microscopic calculation of optical po- 
tentials and proton scattering on 16 [6]. 

An interesting avenue to compute nuclear reactions mi- 
croscopically lies in coupled-cluster theory. The coupled- 
cluster method is a microscopic theory which comes at a 
relatively low computational cost and at the same time 
can provide accurate description of low-lying states and 
properties of nuclei with closed (sub-) shells |8]. The 
coupled-cluster method has in the last decade made sig- 



nificant progress in computing structure of nuclei from 
the valley of stability towards the neutron dripline. 
Using a Berggren basis [9 consisting of bound, reso- 
nant and scattering states, both loosely bound and un- 
bound states have been accurately computed within the 
coupled-cluster formalism [TU1 [TT]. However, so far no 
attempt has been made to apply the coupled-cluster 
method to compute reaction observables, and it is the 
aim of this work to fill this gap and to develop a new 
formalism to compute reaction observables such as elas- 
tic scattering cross sections using micrcoscopic coupled- 
cluster theory. As a first application we will consider the 
elastic scattering reaction 40 Ca(p,p) 40 Ca, whose phase 
shifts and differential elastic cross sections will be evalu- 
ated at low energies. 

Hamiltonian and treatment of the infinite-range 
Coulomb interaction - The intrinsic A— nucleon Hamil- 
tonian consists of kinetic, nuclear and Coulomb parts, 




(1) 

Here, the intrinsic kinetic energy depends on the mass 
number A. The potential Vnn denotes the chiral NN in- 
teraction at next-to-next-to-next-to leading order p~2). [13] 
(with cutoff A = 500 MeV), Vboui is the Coulomb inter- 
action, while V3Ncff is a schematic potential based on the 
in-medium chiral NN interaction by Holt et al. [T3] . The 
potential V^Neff results from integrating one nucleon in 
the leading-order chiral 3NF over the Fermi sphere with 
Fermi momentum kp in symmetric nuclear matter. It 
depends formally on the Fermi momentum kp, the low- 
energy constants cp, and cp of the short-ranged contri- 
butions to the leading-order chiral 3NF, and the chiral 
cutoff. The latter is equal to the value employed in the 
chiral AW interaction [13]. In this work we employ the 
parameters kp = 0.95 £m _1 , cp> = —0.2 and cp = 0.735 
which was recently applied for the study of shell evolution 
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in neutron rich calcium isotopes [T]. 

Let us briefly discuss our treatment of the short and 
long-range parts of the Hamiltonian in Eq. ([IJ. The 
nuclear interaction Vnn is of short range and is ade- 
quately expanded in a basis of harmonic oscillator states 
(see Ref. [TS] for details). The difficulty induced by the 
infinite-range character of the Hamiltonian is thus em- 
bodied in the Coulomb interaction Vboui, asymptotically 
behaving as (Z — l)e 2 /r, with r the distance between the 
isolated proton and the center of charge of the remain- 
ing part of the nucleus. Clearly, it is insufficient to treat 
Vcoui with a harmonic oscillator expansion as we do for 
Vjvjv ■ A solution to this problem has been formulated in 
Ref. [IS]. For this, the Coulomb interaction is rewritten 
as sum of two terms: 



reads: 



V, 



Coul 



U Con \{r) + [Vcoui ~ ^Coui (?")]> 
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where one demands the Coulomb one-body potential 
f^Coui( r ) to behave as [Z — l)e 2 /r for r —> +00. In this 
work we choose £/coui(?0 = erf (or) (Z — l)e 2 /r, where 
erf is the error function and a = it /A fm _ . Thus, the 
[Vboul — f^CouiW] term is short-ranged, so that one can 
use the harmonic oscillator expansion method of Ref. [TS] 
to calculate its matrix elements. Note that the r coor- 
dinate can be taken with respect to the origin of the 
laboratory, because center of charge effects are negligi- 
ble in the asymptotic region on the one hand, and for a 
medium mass nucleus such as 40 Ca on the other hand. 

In order to account for the scattering continuum us- 
ing the coupled- cluster formalism, it is convenient to ex- 
press the Hamiltonian for given partial waves in a basis 
of spherical Bessel functions [11 . Thus, in order to pro- 
ceed, we express Uc OVL \{r) in momentum space, and write 
it in the following way, 



U C oui(k,k') = (k\U Co m(r) - {Z 1)e V ) 



(Z-l)e a _ k 2 + k" 

Qi 



2kk> 
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where £ is the orbital angular momentum of the consid- 
ered partial wave and Qi is the Legendre function of the 
second kind [T7] . As the the first term of Eq. ^ decreases 
very quickly for r — > +00, it can be calculated by numeri- 
cal integration. However, the second term presents a log- 
arithmic singularity at k — kl . In order to counter this 
state of affairs, we will follow the off-diagonal method in- 
troduced in Ref. [TS] ■ It consists in replacing the infinite 
value Qe(l) in Eq. Q occurring at k — k' by a finite 
value depending on the discretization used (see Ref. [18] 
for method and details). 

In order to show the precision of the method in the con- 
text of momentum space calculation, we will diagonalize 
with a basis of Bessel functions the one-body Hamilto- 
nian for the I = partial wave studied in Ref. [18 , which 
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+ EWi(r) (4) 



where m is the proton mass, d = 0.65 fm, Rq — 3 fm, V 
= 52 MeV, and C/couiM is the Coulomb potential of [15]. 
Obtained scattering wave functions have been fitted for 
large r with their asymptotic limit equal to: 



c F(i,r,,kr) + c G(e,r,,kr) 



(5) 



where F(£,rj,x) and G(£,ri,x) are respectively the regu- 
lar and irregular Coulomb wave functions |17j . r\ is the 
Sommerfeld parameter and Cf and Cq are integration 
constants. Regular and irregular Coulomb wave func- 
tions are evaluated numerically using the publicly avail- 
able CWFCOMPLEX code [H| , while C'f and Gq constants 
are determined by fitting Eq. (J5| to the considered scat- 
tering wave functions at r = 10 fm. Results are depicted 
in Fig. [l] It is therein clear that their asymptotic behav- 
ior is very well reproduced, as Coulomb asymptotic ex- 
pansions and diagonalized scattering wave functions are 
virtually indistinguishable for r > 7 fm. This proves that 
the infinite-range character of the Coulomb interaction 
can be handled precisely with Fourier-Bessel transform, 
so that reactions involving protons, such as elastic scat- 
tering, can be undertaken. 
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FIG. 1: (Color online) Scattering s-wave functions <fi(r) ob- 
tained from diagonalization of the one-body Hamiltonian de- 
fined in Eq. Q in momentum space (solid lines) and their 
asymptotic expansion defined in Eq. (JsJ (dashed lines) as a 
function of radius, provided in fm. Wave functions are given 
in units of fm -3 and their energy is written on the figure in 
units of MeV 

One-nucleon overlap functions and coupled- cluster the- 
ory — The scattering of a nucleon on a target A can 
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be described by the one nucleon overlap function. The 
one-nucleon radial overlap function A+1 (lj; r) is defined 
microscopically as the overlap between two independent 
many- nucleon wave functions of A and A + 1 nucleons, 



A+ \lyM) = f (A + l 



A)cp nl] (r). (6) 



The double bar denotes a reduced matrix element, and 
the integral-sum over n represents both the sum over the 
discrete spectrum and an integral over the correspond- 
ing continuum part of the spectrum. The creation op- 
erator aj.. is a spherical tensor of rank j. The radial 
single-particle basis function is given by the term (j) n ij (r) , 
where I and j denote the single-particle orbital and an- 
gular momentum, respectively, and n is the nodal quan- 
tum number. The isospin quantum number has been 
suppressed. The one-nucleon overlap functions describes 
the capture or scattering of an incoming particle with 
quantum numbers Ij on the target nucleus A and with 
the final state A + 1 being either a bound or a scatter- 
ing state. The momentum k is given from the energy 
difference k = y /2rh(E A + 1 - E A )/h, in the case of A 
and A+l being in their ground state E A+1 — E A is 
the one-nucleon separation energy in the A + l nucleus 
(here m = (1 — 1/A)m). We emphasize that the overlap 
function is defined microscopically and independently of 
the single-particle basis. It is uniquely determined by 
the many-body wave functions \ A) and \A + 1). \A) and 
\A+ 1) can in general either be in their ground- or any 
excited state. However, in this work we are interested 
in low-energy elastic scattering, which implies that the 
target nuclues \A) is in its ground state before and af- 
ter the scattering. The one-nucleon overlap functions are 
formally solutions of the Dyson equation, which can be 
written in a Schrodinger like form where the self-energy 
takes the place of a non-local and energy dependent op- 
tical potential [6]. Outside the range R of the optical 
potential, the one-nucleon overlap functions for bound 
A + l states takes (k — in) the form, 



n,i+i/2 

T 
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(7) 



and for A + 1 scattering states (fc > 0), 



0* +L (lj; kr) = B l3 {k) [Fijkr) - tan<Jz(fc)Gz,„(fcr)] . 

(8) 

Here is the Whittaker function, Fi tV and 

Gi tV the regular and irregular Coulomb wave func- 
tions, rj is the Sommerfeld parameter (77 = (Z — 
l)e 2 yJm/2\E\), Cij(in) is the asymptotic normalization 
coefficient (ANC), tan is the I'th partial wave scat- 
tering phase shift at momentum k, and Bij(k) is an ar- 
bitrary normalization constant for the scattering states. 
In order to compute the phase shifts at a given energy, 
it is sufficient to know the one-nucleon overlap function 
A+1 (lj; kr) and it's derivative at a given radius r > R. 



In order to obtain 0^ +1 (ij; kr) we need to solve for the 
ground state of the target nucleus A and the ground- and 
excited scattering states in the residual nucleus A + l. 
The coupled-cluster method is a very efficient tool for 
the computation of ground- and low-lying excited states 
in nuclei with a closed (sub-)shell structure and their 
neighbors. In this work the target nucleus A is a closed 
shell nucleus, and we use the coupled-cluster method 
to compute the ground state of A, i.e. \A) = e T \(j)A)- 
Here is the Hartree-Fock reference state while T is 
a linear combination of particle-hole excitation operators. 
For the residual A + l nucleus we use particle-attached 
equation-of-motion coupled cluster theory to obtain the 
ground- and excited states, and the A + l wave func- 
tions are therefore given by (A + 1| M = ((j) J i L \L A+1 e~ T , 
with L A+1 a linear combination of one-particle, and two- 
particle-one-hole excitations operators (details on our im- 
plementation are presented in Refs. [2U] HI])- Insert- 
ing these expressions for the A and A + l systems into 
Eq. ^, we obtain the coupled-cluster formulation of the 
one-nucleon overlap functions, 



A+l {lj-kr)=t 



nlj 



<f>A) <Pnlj(r), (9) 



here ~ a llj = 



-T~t 



e T is the similarity transformed cre- 
ation operator. The derivation of the diagrammatic and 

algebraic expressions of Eq. (9) and a^y can be found 
in Ref. [22]. Note that in order to compute the radial 
overlap in Eq. ^ we need to use the same mass number 
(41) in the intrinsic kinetic energy of the A and ^4 + 1 
Hamiltonians in Eq. (JIJ. This introduces a small error 
in the ground state of the target nucleus A. However, 
this error decreases rapidly with increasing mass, and we 
estimate that the error is of the order of 100 — 200 keV 
in the relative energy entering the overlap function |23) . 

In our coupled-cluster calculations we use a model 
space consisting of A^ max = 17 major spherical oscilla- 
tor shells with the oscillator frequency hu — 26 MeV. 
This is a sufRccntly large model space to reach practi- 
cally converged results for the ground state of 40 Ca (see 
Ref. pQ). In order to properly account for scattering con- 
tinuum in 41 Sc we use a Gamow-Hartree-Fock basis [24] 
for the relevant proton partial waves. In constructing the 
single-particle basis with the correct treatment of long- 
range Coulomb effects, we use the off-diagonal method in 
momentum space and discretize the one-body momentum 
space Schrodinger equation with 50 mesh points. We find 
that this is a sufficiently large number of mesh points in 
order to obtain the correct Coulomb asymptotics neces- 
sary to describe proton elastic scattering on 40 Ca. 

Results - Figure [2] shows the computed radial over- 
lap function for the ground state of 40 Ca with the J n = 
7/2~ ground state of 41 Sc on a logarithmic scale. Our 
computed proton separation energy for 41 Sc is S^ = 
0.71 MeV, which is in good agreement with the experi- 
mental proton separation energy S^ xp — 1.09 MeV. From 
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the radial overlap function and the separation energy we 
can compute the behaviour of the the overlap function at 
distances beyond the range of the nuclear interaction ac- 
cording to Eq. ([7]). It is clearly seen that the overlap func- 
tion and the known asymptotic form completely overlap 
for distances larger than r ~ 8 fm. Figure [3] shows the 
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FIG. 3: (Color online) Radial overlap functions 0^ +1 (lj; kr) 
between the ground state of 40 Ca and two J 71 = 7/2~ 
scattering states in 41 Sc (solid lines), also shown 
are the corresponding Coulomb scattering functions 
Bij(k) [Fi iV {kr) — t&nSi{k)Gi^{kr)] (dashed lines). 



FIG. 2: (Color online) Radial overlap function 0^ +1 (Ij; kr) 
between the ground state of 40 Ca and the J" = 7/2~ ground 
state of 41 Sc (solid line), also shown is the corresponding 
Whittaker function CijW- v ,i+i/2(iKr)/r for the /r/2 proton 
partial wave (dashed line). 

computed radial overlap functions for the ground state 
of 40 Ca with two — 7/2~ scattering states of 41 Sc, 
at the energies E = 5.439 MeV and E = 16.304 MeV, 
respectively. As we found for the bound overlap function 
shown in Fig. [2j we see that the radial overlap function 
for scattering states and the known asymptotic forms 
completely overlap for distances larger than r ~ 8 fm. 
By matching the asymptotic forms of the overlap func- 
tions given in Eq. Q with the computed overlap func- 
tions, it is clear that we can determine the correspond- 
ing elastic scattering phase-shift at the computed scat- 
tering energy. Figure [4] shows our computed scattering 
phase shifts for proton elastic scattering on 40 Ca for the 
s%/2,Pi/2,P3/2, ^3/2)^5/2 partial waves at energies below 
14 MeV. The solid dots correspond to the computed scat- 
tering energies, and we used cubic spline to interpolate 
between the discrete set of scattering energies. We clearly 
see the appearance of a narrow resonance in the p 3 / 2 par- 
tial wave around ~ 1.6 MeV, while a broader resonance 
appear in the pi/2 partial wave at around ~ 3.4 MeV. 
In order to check our results we computed the low-lying 
resonances in 41 Sc using a complex Gamow-Hartree-Fock 
basis [IT], and we found a J" = 3/2~ resonance at the 
energy E = 1.61 - O.OOli MeV and a ,V = 1/2" res- 
onance at the energy E = 3.42 - 0.20i MeV. Clearly 



these energies are consistent with the resonances appear- 
ing in the P3/2 and P1/2 elastic scattering phase-shifts 
in Fig. [4] From the scattering phase-shifts we can com- 
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FIG. 4: (Color online) Computed phase-shifts for elastic pro- 
ton scattering on 40 Ca for low-lying partial waves and energies 
below 14 MeV. 

pute the differential cross section for elastic proton scat- 
tering as described in e.g. Ref. 25j. Figures [5] and § 
show the differential cross section divided by the Ruther- 
ford cross section for elastic proton scattering on 40 Ca 
at the relative-center-of mass energies E cm = 9.6 MeV 
and E cm = 12.44 MeV, respectively. All partial waves 
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for I < 2, were included in the computation of the cross 
sections. Overall we get good agreement between our 
calculated cross sections and the experimental cross sec- 
tions. In particular we sec that our computed minima are 
in good agreement with the experimental minima, while 
we tend to overestimate the cross sections at large scat- 
tering angles. The overestimated cross sections at large 
angles is most likely due to the fact that we do not ac- 
count for intermediate excitations that takes place above 
the deuteron threshold, and these excitations generally 
cause absorption and reduces the cross section at large 
angles. Going beyond 2p-lh excitations for the compu- 
tation of the A + 1 wave functions will account for such 
effects, and we are working towards such improvements 
in our approach. We also computed the cross sections 
including the /s/2, It/2-, 37/2; 59/2 partial waves, however 
the agreement with data did not improve. This is most 
likely due to the fact that our computations for the f§/2 
and gg/2 partial waves finds very narrow resonances at 
too low energy energies as compared to experiment. This 
is also consistent with the ~ 10 MeV overbinding we get 
for 40 Ca using the Hamiltonian in Eq. (fiT) pQ. 
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FIG. 6: (Color online) Differential cross section divided by 
Rutherford cross section for elastic proton scattering on 40 Ca 
at -E cm = 12.44 MeV(solid line), experimental data (dots) is 
taken from Ref . [26] . 
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FIG. 5: (Color online) Differential cross section divided by 
Rutherford cross section for elastic proton scattering on 40 Ca 
at -Ecm = 9.6 MeV (solid line), experimental data (dots) are 
taken from Ref. [26] . 

Conclusions — Using coupled-cluster theory, we com- 
puted cross sections for elastic scattering of protons 
on 40 Ca, at the center-of-mass energies 9.6 MeV and 
12.44 MeV, respectively. We found good agreement for 
our computed diffraction minima with experiment, while 
we tend to overestimate the cross sections at large scat- 
tering angles. The key ingredients for computing observ- 
ables for proton scattering are; (i) the one-nucleon over- 
lap function computed from microscopic coupled-cluster 
theory, and (ii) a single-particle basis that has the correct 



Coulomb asymptotics. We showed that the newly devel- 
oped off-diagonal method is a very accurate method to 
compute Coulomb scattering wave functions in momen- 
tum space. The fast convergence of the scattering wave 
functions with increasing number of mesh-points makes 
this basis an ideal starting point for computing reaction 
observables. This work constitute the first successful ap- 
plication of coupled-cluster theory to nuclear reactions, 
and we believe it makes a significant leap forward in link- 
ing reactions with microscopic structure calculations. 
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